Insect Biochemistry and Molecular Biology

The European honey bee (Apis mellifera L.) is a key agricultural pollinator frequently exposed to pesticide residues, yet the molecular basis of its chemical adaptation, particularly glutathione S-transferases (GSTs) involved in xenobiotic detoxification, remain incompletely understood. In this study, AmGSTO1 was structurally and functionally characterized to evaluate its role in agrochemical interaction and protection against oxidative stress. The crystal structure of AmGSTO1 in complex with glutathione revealed its 3D architecture and key active-site residues were identified by structural analysis and site-directed mutagenesis. Fluorescence binding assays demonstrated measurable affinity for multiple agrochemicals, including TCP, fenoprop, 2,4-D, tetramethrin, nicotine, and 3-phenoxybenzaldehyde. However, HPLC analysis showed no detectable substrate depletion, suggesting ligand binding to AmGSTO1 without catalytic turnover. AmGSTO1 exhibited antioxidant activity toward cumene hydroperoxide, hydrogen peroxide, and paraquat, as well as dehydroascorbate reductase activity. These findings indicate that AmGSTO1 may contribute to agrochemical tolerance through ligand sequestration and redox protection mechanisms.

Adult females of the two-spotted spider mite, Tetranychus urticae Koch, enter a photoperiodically induced diapause to overwinter. Diapause in T. urticae is accompanied by reproductive arrest and the orange body coloration that arises from the accumulation of astaxanthin esters. How these two traits are coordinated at the molecular level remains poorly understood. Here, we compared the proteomes of adult females reared under diapause-inducing (long-night) and non-diapause-inducing (short-night) photoperiods using liquid chromatography-tandem mass spectrometry, followed by RNA interference (RNAi) of candidate genes. The carotenoid biosynthesis enzymes phytoene desaturase (TuPDS) and lycopene cyclase/phytoene synthase (TuLCPS), both encoded by genes horizontally transferred from fungi, were more abundant in diapausing females than in non-diapausing females. RNAi of the genes encoding TuPDS and TuLCPS markedly reduced orange pigmentation as well as {beta}-carotene and astaxanthin contents, demonstrating that these enzymes are required for diapause-associated pigmentation. Our proteomic analysis further identified a single PLAT (Polycystin-1, Lipoxygenase, Alpha-toxin) domain protein, TuPLAT10, as one of the most strongly upregulated proteins in diapausing females. The PLAT domain is a lipid-binding module, suggesting a role for TuPLAT10 in lipid metabolism. In addition to the suppression of orange pigmentation, RNAi of the TuPLAT10 gene restored reproduction even under diapause-inducing conditions and selectively reduced TuPDS and TuLCPS protein levels, despite the absence of sequence similarity to their genes. We propose that TuPLAT10 acts as a lipid-allocation switch that, in response to photoperiodic information, partitions fatty acids between astaxanthin esterification and yolk lipid supply, thereby coupling reproductive arrest and carotenoid pigmentation during diapause in T. urticae.

A proteomic analysis of Ixodes scapularis nymph saliva identified 252 proteins, including six tubular lipid-binding proteins (TULIPs). Comparing nymphs fed on mice that were uninfected or infected with Borrelia burgdorferi, twelve salivary proteins showed significant differences in the amounts detected, including XP_040079658.2, which we refer to as TULIP2. Considering the known immunity-related functions of some TULIPs, we expressed and purified TULIP2 from Escherichia coli and analyzed its interaction with B. burgdorferi lipids. The purification of TULIP2 from E. coli presented many obstacles, due to insolubility, which is consistent with previous reports from studies of other TULIP family members. The binding results showed specificity for B. burgdorferi lipids, with evidence for cholesteryl {beta}-galactoside as a major binding target. Molecular modeling of TULIP2 did not show any strong lipid binding sites. We used molecular dynamics simulation of TULIP2 to explore its conformational landscape by thermal unfolding. The earliest unfolding intermediate opened a hydrophobic pocket to which cholesteryl {beta}-galactoside was predicted to bind strongly. We propose that a specific lipid bilayer interaction with TULIP2 triggers the opening of the ligand-binding site.

CRISPR-Cas genome editing toolkits have expanded the scope of genetic studies in various emerging model organisms. However, their applications are limited mainly to knockout experiments due to technical difficulties in establishing knock-in strains, which enable in vivo molecular tagging-based experiments. Here, we investigated knock-in strategies in the harlequin ladybug Harmonia axyridis, a model insect for evolutionary developmental biology, which shows more than 200 color pattern variations within a species. We tested several knock-in strategies using synthetic DNA templates. We found that ssDNA templates generated founder knock-in strains efficiently (2.5-11%), whereas the 5 regions of ssDNA templates were frequently deleted when the insert length exceeded [~]40 bases. To overcome this limitation, we designed several 3 extended DNA templates. Fast-annealed 3-extended double-stranded DNA templates, which were designed for tagging endogenous proteins with epitope tags, showed high founder generation efficiency (9.9-20.9%) and accuracy (30.8-85.7%). This strategy is also applicable to the two-spotted cricket Gryllus bimaculatus, suggesting that the fast-annealed 3-extended dsDNA template is a versatile DNA template for generating knock-in strains in emerging model insects for developmental genetic studies. Summary statementFast-annealed 3-extended dsDNA templates facilitate efficient CRISPR-Cas9-mediated knock-in in emerging model insects.

The brain is one of the most complex animal organs but the development of the many different neuron types remains enigmatic. A set of brain-specific transcription factors is known to be involved in brain patterning but their specific contributions remain to be elucidated in most cases, including foxQ2II. This transcription factor is known to be conserved in anterior neuroectodermal patterning of most animals while it has been lost from vertebrates. However, the contribution of foxQ2II-positive neurons to the adult brain has remained enigmatic. Here, we use an enhancer trap, immunostainings and our newly established beetle brainbow system to categorize Tc-foxQ2II-positive neurons into nine clusters with different projection patterns. All clusters contain neurons with the fast activating neurotransmitters acetylcholine and glutamate while no Tc-foxQ2II positive neuron is GABA-ergic or serotonin-positive. Interestingly, we found that many dopaminergic neurons were Tc-foxQ2II positive and we homologize them with dopaminergic neurons of the PPL2c, PPM1 and PPL1 cluster described in the Drosophila brain. Our results show that Tc-foxQ2II marks subsets of fast-acting interneurons contributing to the higher order brain centers mushroom bodies and central complex. Taken together, our work expands the known functional range of foxQ2 genes from sensory and neurosecretory cell specification to interneurons involved in the function of higher order brain centers.

Octopamine is involved in a variety of different physiological and behavioral mecha-nisms in Drosophila melanogaster. Throughout the life cycle of the fruit fly, from the larva to the adult, octopaminergic neurons in both the central and the peripheral nerv-ous system target a multitude of neurons and even non-neuronal tissues, making it challenging to analyze individual mechanisms of octopamine function. One approach to deconstructing this complex system is to examine the postsynaptic components of signal transmission. In Drosophila, octopamine interacts with six distinct G-protein-coupled receptors. For some of these receptors, expression maps and functional im-plications have been described. In contrast, other receptors have been neglected, partly due to the lack of suitable genetic tools. Here, for the first time, we compiled a complete set of mutant lines of all known octopamine receptors, all generated using the same genetic tool, the recently established Trojan Exon system. It integrates the Gal4/UAS binary expression strategy while simultaneously impairing receptor func-tion. This enabled us to generate a comprehensive anatomical map of receptor ex-pression in the larva and, at the same time, analyze the function of individual octopa-mine receptors during larval development, chemosensory perception and locomotion. All octopamine receptors (Oamb, Oct2R, Oct{beta}1R, Oct{beta}2R, Oct{beta}3R, and Oct-TyrR) showed extensive signal in the central nervous system. The same was found for the peripheral nervous system, with the exception of Oct{beta}2R, which showed pronounced expression in the somatic muscles. We also observed a previously undescribed role of Oct{beta}1R, Oct{beta}3R, and Oct-TyrR in larval hatching and in the survival of larvae and pupae. Molecular evaluation of the Trojan Exon octopamine lines supports our analy-sis. In addition, we combined the experimental results with gene expression data from the different development stages of Drosophila melanogaster and from different tis-sues and cell populations throughout the body. Overall, we compiled, analyzed and validated a complete set of octopamine lines which, together with gene expression analysis, provides a basis for further functional studies on the larval octopaminergic system.

Background and AimsPlants deploy triterpenoid saponins as chemical defences against herbivores, yet it remains unclear whether insect digestion detoxifies these compounds or generates equally or more active metabolites. Because saponin bioactivity depends strongly on glycosylation patterns, we examined the fate and defensive activity of hederagenin-derived saponins during herbivory. MethodsLarvae of Plutella xylostella were fed leaf discs containing structurally defined hederagenin-derived saponins. Saponin composition in treated leaves and larval frass was analysed by LC- qTOF-ESI-MS/MS. Feeding assays were used to compare the antifeedant activity of mono- and bidesmosidic forms. Key ResultsLarvae selectively metabolized complex hederagenin-derived saponins into simpler forms, with cellobiosides converted into monoglucosides during digestion, resulting in a marked shift in saponin composition between ingested material and frass. Feeding assays showed that monodesmosidic saponins strongly deterrer feeding, whereas bidesmosidic saponins were largely inactive. The loss of activity in bidesmosidic saponins was not explained by differential metabolism, indicating that glycosylation patterns directly determine biological function. ConclusionsInsect herbivores selectively modify saponin structures through deglycosylation, thereby altering their defensive properties. Our findings demonstrate that glycosylation governs both saponin activity and metabolic fate, highlighting insect-driven turnover as a critical component of plant chemical defence during plant-herbivore interactions. Issue SectionOriginal article

Animals sense and integrate complex external cues to make developmental decisions that help them better survive and adapt to their natural habitats. Under environmental adversity, nematodes can enter an alternative developmental pathway to form a diapautic and stress-resistant stage, termed the dauer larvae. While dauer formation has been well characterized in Caenorhabditis elegans, how environmental factors influence analogous stages in other nematode species remains largely unexplored. This study examines how symbiotic bacteria, temperature, and pheromones affect the formation of the infective juvenile (IJ), a dauer-like stage, of the insect-parasitic nematode Steinernema hermaphroditum. In contrast to C. elegans, where dauer entry is promoted by heat, IJ development in S. hermaphroditum development is enhanced by reduced temperature. Moreover, the presence and absence of live symbiotic bacterium Xenorhabdus griffiniae functions as an ON-and-OFF switch that regulates the host IJ formation. Crude pheromone extracts from S. hermaphroditum liquid culture do not robustly induce IJ formation in a dose-responsive manner, unlike the potent pheromone-driven dauer entry observed in C. elegans. Nutrient-rich liver-kidney media that mimics host insect environment showed IJ entry induction in a pheromone-dependent manner. These data suggest that external cues, such as temperature, microbial diet, and pheromone, are perceived differently by S. hermaphroditum in comparison to that of C. elegans, reflecting species-specific adaptations to distinct ecological niches and life history strategies.

The tropical house cricket Gryllodes sigillatus is a major species used in the edible insect farming industry. Despite the rapid expansion of this sector, diagnostic tools for detecting infections in these species remain limited. The lack of validated reference genes compromises the reliability of RT-qPCR-based gene expression analyses, which are essential for the development of molecular tools for disease diagnosis and health monitoring in insect production systems. To address this gap, we evaluated the expression stability of six candidate reference genes (ACTB, EF1, GAPDH, HisH3, RPL5, and 18SrRNA) across four body parts (abdomen, head, legs, and whole body) using a combination of complementary statistical approaches, including geNorm, NormFinder, BestKeeper, the {Delta}Ct method, the R statistical environment, and the integrated RefFinder tool. Candidate genes were identified and annotated using the recently published G. sigillatus genome, through sequence comparisons with closely related insect species using BLAST and reciprocal BLAST analyses, multiple sequence alignments. All procedures complied with MIQE 2.0 guidelines to ensure methodological rigor and transparency. The results showed that ACTB, EF1, RPL5, and 18SrRNA exhibited stable and consistent expression across all analyzed tissues, whereas GAPDH and HisH3 displayed high variability and were generally unsuitable for normalization, except in head tissue where GAPDH remained stable. This study provides the first validated set of reference genes for G. sigillatus, establishing a robust foundation for accurate, reproducible, and comparable gene expression analyses. Furthermore, these findings support the development of RT-qPCR-based diagnostic tools, contributing to improved health monitoring and biosafety in insect production systems.

The western honey bee (Apis mellifera) is an essential pollinator facing unprecedented threats from pesticide exposure. While pesticide resistance evolution is well documented in agricultural pests, our understanding of genetic variation in honey bee detoxification systems remains limited. This represents a missed opportunity, as harnessing naturally occurring detoxification diversity could provide new avenues for pollinator protection. Cytochrome P450 monooxygenases (CYPs), which are central to xenobiotic metabolism, offer a promising starting point. Here, we present the first comprehensive analysis of CYP genetic diversity in A. mellifera. We analysed the CYPome of 1,467 individuals representing 18 A. mellifera subspecies from 25 countries and identified 5,756 single-nucleotide polymorphisms (SNPs) in 46 CYP genes. Imputed McDonald-Kreitman testing revealed that 56% of non-synonymous CYP substitutions were driven by positive selection. Of the 1,302 haplotypes identified, 84% resided in CYP3, concentrated in the CYP9 and CYP6AS subfamilies implicated in xenobiotic detoxification. Population-level analysis of nucleotide diversity, Tajimas D selection signatures, FST-based differentiation, and McDonald-Kreitman testing pointed to CYP3 clan genes as the primary locus of adaptive variation. This work provides the first step toward building a comprehensive pharmacogenomic resource for honey bees, enabling the prediction of population-specific pesticide vulnerabilities and leveraging naturally occurring detoxification variants to enhance pollinator resilience - a critical step toward sustainable pollinator management.

Nociception is an essential response for organisms to avoid potential harm and promote survival. Its molecular determinants are largely conserved across Eumetazoa. TRPV receptors are polymodal ion channels exhibiting selective peripheral expression and functional coupling that underpins nociception and pain modulation in complex organisms. However, the execution of protective behaviours triggered by TRPVs is also found in species with a simpler nervous organisation, thus encouraging their investigation in invertebrate model organisms to increase understanding of animal nociception. Cephalopods represent an interesting invertebrate phylum with respect to the evolution of the nervous system, whose complexity suggests it might support pain-like states that exist in vertebrates. This possibility is reflected by the inclusion of cephalopods in the UK and EU animal welfare legislations. Despite this, there is poor characterisation of cephalopod molecular nociceptors. For this reason, we used in silico analysis to identify two TRPV channels in Octopus vulgaris genome (Ovtrpv1 and Ovtrpv2). We validated the putative transcript sequences and highlighted prevalent expression in sensory tissues. We investigated the functional competence of these TRPVs by heterologously expressing Ovtrpv1 and Ovtrpv2 cDNA into Caenorhabditis elegans null mutants of the orthologous genes, ocr-2 and osm-9 respectively. Ovtrpvs successfully rescued the aversive response to chemical and mechanical noxious stimuli in the C. elegans mutants, suggesting these receptors are polymodal nociceptors. Additionally, complementary investigation using Xenopus laevis oocytes showed Ovtrpv1 and Ovtrpv2 form an active heteromeric channel gated by nicotinamide. This study highlights Ovtrpvs as an important route to better understand nociceptive detection in cephalopods.

Rhodnius prolixus is a primary insect vector of Trypanosoma cruzi, the causative agent of Chagas disease, a neglected parasitosis endemic in Latin American countries. It has been estimated that Chagas disease affects 7-8 million people worldwide and is responsible for approximately 1000 deaths per year. Genetic and molecular studies in this species remain challenging due to its life cycle and feeding habits, thus hindering the development of new strategies to control their populations and reduce the diffusion of Chagas disease. Recently, two stable cell lines - RPE/LULS53 and RPE/LULS57 - were derived from Rhodnius embryos, which represent promising new tools to investigate the genetics of this insect vector. Here, we describe their gene expression landscapes through transcriptomic approaches. We show that 8,968 expressed genes are shared between the two cell lines, whereas 391 and 1,088 genes are uniquely expressed in RPE/LULS53 and RPE/LULS57, respectively. Although key components of primary developmental, immune and redox signaling pathways are expressed in both cell lines, some genes such as Frizzled-10-a-like and catalase show marked differences in expression. Our results strongly suggest that RPE/LULS53 and RPE/LULS57 likely represent two different cell phenotypes. Consistent with this, gene ontology analysis reveals that RPE/LULS53 is enriched for animal organ morphogenesis and stress response, while RPE/LULS57 for DNA-directed RNA polymerase activity, among others. Despite these differences, both cell lines express comparable levels of transcripts from resident transposable elements, including the highly abundant Mariner and LINE/I elements, as well as horizontally transferred transposons. Our findings shed light on the nature of the RPE/LULS53 and RPE/LULS57 embryo-derived cell lines and provide valuable transcriptomic resources for future genetic and functional studies in Rhodnius and other triatomine insect vectors. Author summaryRhodnius prolixus is a blood-feeding insect and a major vector of Chagas disease, a parasitosis endemic in Latin America and affecting millions of people worldwide. In the absence of effective drugs and vaccines, the control of the insect population represents a promising strategy to reduce the diffusion of the disease. Yet, genetic and functional studies in Rhodnius are extremely challenging due to its feeding habit and life cycle. To overcome these limitations, researchers have previously developed two stable cell lines derived from Rhodnius embryos. In this study, we provide the first characterization of the genes expressed in these cell lines. We found that, while the two cell lines share many expressed genes, each of them also has distinct gene expression patterns pointing to two different cell types with specialized functions. These differences likely affect the way they respond to stress and regulate biological processes. Our findings provide an important resource for researchers studying Rhodnius prolixus and other insect vectors, helping advance our understanding of the genetic and molecular mechanisms that control the insect development and mediate the interactions between insect vectors and the parasites they transmit

Insecticide resistance is threatening malaria control. While the evolution and spread of resistance has been linked to scale-up in the distribution of public health insecticides, the role of environmental pollutants such as the polyaromatic hydrocarbons (PAHs) from industrial and agricultural use remains largely uncharacterized. The PAHs are potent ligands of the aryl hydrocarbon receptor (Ahr) transcription factors involved in the regulation of xenobiotic metabolizing enzymes, and potentially involved in insecticide resistance. Here, using field insecticide-resistant (Auyo) An. coluzzii and a laboratory-susceptible colony (Ngousso), we conducted a multi-generational selection experiment using naphthalene, fluorene and a mixture of both PAHs. After ten generations, the changes in susceptibility to insecticides were monitored using WHO bioassays and whole-transcriptome analysis (RNASeq) was conducted. Compared with the non-selected colony lines, PAH exposures significantly reduced pyrethroid and DDT resistance in the field population, suggesting fitness cost associated with established resistance. In contrast, Ngousso showed a significant increase in DDT resistance (p = 0.01) at the tenth generation. A significant increase in permethrin resistance was also observed at the seventh generation (p = 0.03). Several candidate genes from the major detoxification classes were overexpressed in the selected lines (including GSTe2, CYP6Z1, and CYP6P4); the most consistent were CYP6M4 and CYP4C27, as well as those from the Ahr pathway. Heterologous expression of CYP6M4 revealed its ability to metabolise pyrethroids, including permethrin, deltamethrin, and -cypermethrin, as well as PAHs (naphthalene and fluorene). These findings establish the role of environmental pollutants as additional drivers of metabolic insecticide resistance in An, coluzzii.

The black legged tick, Ixodes scapularis, is a vector of the bacterium that causes Lyme disease and several other illnesses, including anaplasmosis, babesiosis, and tick-borne encephalitis. Although high-quality genome annotations are available for I. scapularis, functional understanding of I. scapularis genes is limited. To address this, we developed a platform for genome-wide CRISPR-Cas9 knockout screening in I. scapularis cells. To evaluate the platform, we performed a screen to identify genes associated with cellular fitness, and screens for resistance to treatment with copper chloride, Antimycin A, or Destruxin A (DA), a cyclic hexadepsipeptide produced by the pathogenic fungus Metarhizium anisopliae. In each case, the screens implicate specific sets of conserved and non-conserved I. scapularis genes in relevant cellular functions, providing the first experimental evidence of function for a large set of I. scapularis genes. Altogether, in this first-of-its-kind effort for the arthropod subclass Acari, we present an unbiased genome-wide CRISPR-Cas9 knockout cell screening platform, related resources, and datasets that will be broadly useful to efficiently uncover cellular functions of I. scapularis genes.

The cockroach Blattella germanica possesses panoistic ovaries, in which oocytes lack nurse cells and therefore need to rely on their own transcriptional activity to support embryogenesis. Ovarian development in this species involves the development of a single basal ovarian follicle (BOF) per gonadotropic cycle, a process strictly regulated by endocrine signals, primarily juvenile hormone and ecdysone, which act at both the transcriptional and translational levels. In addition, transcriptional activity in these ovaries is necessary for both regulating and genome protection, and at this level, PIWI-interacting RNAs (piRNAs) play an essential role. Although insect ovaries are known to be particularly rich in piRNAs, their function in ovary maturation is still not well defined. For this purpose, we characterize the piRNA expression dynamics across seven key developmental and reproductive stages, ranging from late nymphal instars to post-vitellogenic adults. piRNA expression in B. germanica shows coordinated fluctuations. Expression remains stable in previtellogenic ovaries, whereas vitellogenic ovaries show pronounced changes. Moreover, vitellogenic ovaries exhibit reduced piRNA diversity due to strong enrichment of a subset of highly expressed piRNAs. Our data show that although piRNAs predominantly map to transposable elements, particularly LINEs, there is a notable increase in gene-derived piRNAs toward the end of the cycle. Our results suggest regulatory roles of piRNAs in modulating both TEs and mRNAs during BOF maturation, likely related to changes in the follicular cell program.

To improve their chances of reproductive success, nematodes not only must arrest their development in response to adverse growth conditions but also must quickly recover if conditions improve. A polyketide synthase (PKS)-nonribosomal peptide synthetase (NRPS) hybrid assembly line that is expressed in the canal-associated neurons (CANs) of Caenorhabditis elegans promotes recovery from starvation-induced larval arrest. Here, we show that in the predatory nematode Pristionchus pacificus this assembly line produces a suite of secondary metabolites, including a family of hybrid polyketide-nonribosomal peptides known as the nemamides, the related nematides, and a family of ascarylose-modified polyketides named ascarenes. Depending on the starter unit that is loaded onto the PKS, the assembly line can produce dramatically different downstream products. Whereas the nemamides promote recovery from starvation-induced larval arrest, the ascarenes inhibit development of the dauer larval stage and promote recovery. This dichotomy suggests that the PKS-NRPS megasynthetase serves as a signaling hub in the CANs, controlling multiple developmental events. The PKS-NRPS assembly line is highly conserved across many nematode species, and identification of these chemical signals will help to elucidate the signaling pathways that control development in the worm and lead to novel anthelmintics.

MicroRNAs (miRNAs) play essential roles in the posttranscriptional regulation of gene expression in organisms. In the process of synthesizing mature miRNAs from miRNA precursors, the miRNA precursors are cleaved via Dicer at their loop structure, after which the miRNA precursors become mature and regulate transcription. However, the consequences of altering the loop sequence are not fully understood. The silkworm Bombyx mori is a lepidopteran insect with many genetic strains. We identified a mutant of the miRNA miR-3260 whose the part of the loop structure was lacking in a silkworm strain with translucent larval skin. Here, we aimed to analyze the role of wild-type miR-3260 and the influence of the mutation of the loop structure in B. mori. First, we identified the genomic region responsible for the translucent larval skin phenotype and determined that the mutated miR-3260 nucleotide sequences. Then, we predicted the binding partners of wild-type miR-3260 using the RNA hybrid tool and found two juvenile hormone (JH)-related genes as targets of wild-type miR-3260. Next, we assessed the relationships between miR-3260 and JH and found that miR-3260 was highly expressed in the Corpora allata and its expression responded to JH treatment. Meanwhile, miR-3260 mimic and inhibitor did not induce the typical phenotypes associated with JH in B. mori. Then, we compared the dicing products from wild-type and mutant miR-3260 precursors and observed that neither form underwent Dicer-mediated cleavage when the loop structure was altered. These results suggest that loop mutations in the miR-3260 precursor may not influence dicing activity, consistent with the lack of observable phenotypic effects.

Gene model for the ortholog of Lst8 (Lst8) in the May 2011 (WUGSC dyak_caf1/DyakCAF1) Genome Assembly (GenBank Accession: GCA_000005975.1) of Drosophila yakuba. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

BACKGROUNDThe two-spotted spider mite, Tetranychus urticae Koch, is a major agricultural pest with a rapid propensity for developing acaricide resistance. Bifenazate targets mitochondrial cytochrome b (CYTB). While the G126S mutation is frequently associated with resistance, its independent role remains unclear as it often occurs with other substitutions. This study explores the molecular basis of bifenazate resistance in a Russian laboratory strain derived from a St. Petersburg greenhouse population. RESULTSDisruptive selection with increasing bifenazate concentrations generated resistant and susceptible isofemale lines. AlphaFold2 structural modeling of CYTB indicated that G126S causes a steric clash, leading to conformational destabilization, whereas other reported mutations primarily affect the ligand-binding pocket. Oxford Nanopore sequencing revealed a very low initial frequency of the G126S allele (<1%; 226/35,895 reads) in the unselected population. After one year of stepwise selection (0.00005-0.031% a.i.), the mutant allele frequency surged to 90% (7,272/8,056 reads). No other known resistance-associated mutations were found in the analyzed cytb fragment. CONCLUSIONWe report the first identification of the G126S mutation in a Russian T. urticae population and demonstrate its rapid fixation under bifenazate selection. Within this genetic background, G126S alone appears sufficient to confer high-level resistance, emphasizing the population-specific nature of resistance evolution and the critical need for local monitoring.

The Drosophila adult central brain contains 240 circadian neurons, of which there are more than 25 different neuron subtypes based on connectomic data. Recent single cell RNA-seq (scRNAseq) characterization of these neurons "around the clock" also indicates a similar number of molecular subtypes of circadian neurons, but other conclusions from these transcriptomic studies warranted verifying and extending with other approaches. To this end: 1) We used a genetic multiplexing strategy to profile the transcriptomes of circadian neurons from multiple time points in a single experiment, reducing confounding technical variation between timepoints; 2) Large numbers of single nuclei were sequenced (snRNA-seq), which was enabled because the new method EL-INTACT purifies nuclei from frozen heads; 3) We assayed 12 time points under both light-dark (LD) and constant darkness (DD) conditions. These approaches showed dramatic transcriptional differences between time points in many circadian neuron types and enhanced time-of-day gene expression analysis. The data indicate that most of this regulation is transcriptional and circadian. There were however a small number of light-dependent transcripts, including a few that correspond to mammalian immediate-early genes. They probably play a role in the light-regulation of gene expression and behavior in specific neurons, perhaps circadian entrainment or phase-shifting. The results taken together provide a more comprehensive picture of gene expression heterogeneity within adult Drosophila circadian neurons including how intrinsic clock mechanisms and light cues are integrated across circadian neuron subtypes.